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Self-induced flavor conversions of supernova (SN) neutrinos have been characterized in the spher¬ 
ically symmetric “bulb” model, reducing the neutrino evolution to a one dimensional problem along 
a radial direction. We lift this assumption, presenting a two-dimensional toy-model where neutrino 
beams are launched in many different directions from a ring. We find that self-interacting neutrinos 
spontaneously break the spatial symmetries of this model. As a result the flavor content and the 
lepton number of the neutrino gas would acquire seizable direction-dependent variations, breaking 
the coherent behavior found in the symmetric case. This finding would suggest that the previous 
results of the self-induced flavor evolution obtained in one-dimensional models should be critically 
re-examined. 

PACS numbers: 14.60.Pq, 97.60.Bw 


I. INTRODUCTION 

Dense neutrino gases in early universe or emitted from 
core-collapse supernovae (SNe) represent unique cases to 
probe the effect of the neutrino-neutrino interactions on 
the flavor conversions. Indeed, in these environments 
the neutrino-neutrino interactions would generate a large 
neutrino potential ^ ~ y/^Cpni, that in some cases can 
exceed the ordinary matter term A = y/2Gpne and the 
neutrino vacuum oscillation frequency w = Am'^/2E. 
When this situation is encountered the neutrino-neutrino 
potential would dominate the flavor evolution producing 
large self-induced flavor conversions (see [l| for a review). 
A vivid activity on these effects in the context of SN neu¬ 
trinos has flourished since a decade Indeed, it has 

been realized that in the deepest SN regions self-induced 
effects can produce collective neutrino oscillations, lead¬ 
ing to peculiar spectral features in the oscillated neutrino 
spectra, dubbed as spectral swaps and splits US]- 

The development of the self-induced flavor conversions 
is associated with instabilities in the flavor space that are 
triggered by the interacting neutrinos. The first one to 
be noticed was the himodal instability present even in an 
homogeneous and isotropic neutrino gas Q. In particu¬ 
lar, it was shown that an ensemble initially composed of 
equal densities of Ve and De in the presence of a dominant 
neutrino-neutrino interaction term would exhibit in in¬ 
verted mass hierarchy (Am^ < 0) large pair-conversions 
of the type VeUe ^ even with a small mixing angle. 
This behavior has been explained in terms of an unstable 
pendulum in flavor space, where the instability is associ¬ 
ated with the tiny mixing angle aini- Furthermore, it 
has been shown that if one introduces an anisotropy in 
the neutrino gas, this can dramatically change the pre¬ 
vious solution. Indeed, in a non-isotropic neutrino en¬ 
semble, the neutrino-neutrino interaction term contains 
multi-angle effects since the current-current nature of the 
low-energy weak interactions introduces an angle depen¬ 
dent term (I — Vp • Vq) between two interacting neutrino 


modes i[ll. In the case of a gas completely symmet¬ 
ric in flavor content of v and P, even a small deviation 
from a perfect isotropy is enough to produce a multi¬ 
angle dechoerence leading to a flavor equilibrium among 
the different neutrino species in both the mass hierar¬ 
chies M- Multi-angle effects have been extensively stud¬ 
ied in the context of flavor evolution of SN neutrinos M, 
whose emission is far from isotropic. It has been realized 
that in some cases they can destroy the collective behav¬ 
ior of the flavor evolution observed in an isotropic envi¬ 
ronment Multi-angle effects can also lead 

to a trajectory-dependent matter term, which if strong 
enough suppresses the self-induced conversions [T7l - [20l |. 
In the context of SN neutrinos it has been often assumed 
an axially symmetric neutrino emission in oder to inte¬ 
grate out the azimuthal angle in the multi-angle kernel. 
However, it has been found that lifting this assumption, 
neutrino-neutrino interactions can break axial symme¬ 
try and lead to azimuthal-angle dependent flavor conver¬ 
sions [2 iMT7I| . 

The lesson that has been gained from these situations 
is that self-interacting neutrinos can spontaneously break 
the symmetries of the initial conditions, since small de¬ 
viations from them can be dramatically amplified during 
the further flavor evolution. This insight has recently 
stimulated doubts about the validity of the solution of 
the SN neutrino equations of motion worked out in the 
so-called “bulb model” [ 1 , [l^. In this framework it is 

assumed the spherical symmetry about the center of the 
SN and the axial symmetry about any radial direction. 
These two symmetries allow one to reduce the problem 
to a one-dimensional evolution along a radial direction. 
Remarkably, removing the assumption of spherical sym¬ 
metry it would necessary to solve a challenging multi¬ 
dimensional problem to characterize the neutrino flavor 
evolution. 

In this context, in order to show how deviations from 
the spatial symmetries of a system would affect the fla¬ 
vor evolution a simple two-dimensional model has been 
recently proposed in Namely, monochromatic neu- 
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trinos streaming in a stationary way in two directions 
(“left” L and “right” R, respectively) from an infinite 
boundary plane at z = 0 with periodic conditions on 
X and translation invariance along the y direction. Re¬ 
markably, there is a correspondence between the symme¬ 
tries of the bulb model and the ones of this planar case. 
Indeed, the translational symmetry in the x direction in 
the planar model corresponds to the spherical symmetry 
of the bulb-model and the L-R symmetry is equivalent 
to the axial symmetry in the spherical case. By means 
of a stability analysis of the linearized equations of mo¬ 
tion, it has been shown in the planar model that if one 
perturbs the initial symmetries of the flavor content in 
both the two emission modes and along the boundary in 
the X direction, then self-induced oscillations can sponta¬ 
neously break both these spatial symmetries . In [2^ 
we have recently performed a numerical study of the fla¬ 
vor evolution for this case. We found that the initial 
small perturbations are amplified by neutrino interac¬ 
tions, leading to non-trivial two-dimensional structures 
in the flavor content and lepton number of the neutrino 
enseble, that would exhibit large space fluctuations. 

The purpose of this paper is to develop a two- 
dimensional model to capture more closely some of the 
features of the SN environment. In particular, with re¬ 
spect to the planar model considered in [H, [2^ we make 
the following improvements: (i) h' emission from a ring 
mimicking the neutrino-sphere, (ii) parameters inspired 
by the SN neutrino emissivity, (in) declining neutrino 
density from the boundary, (iv) multi-angle effects. We 
also assume that self-induced flavor conversions would 
develop without any hindrance caused by a large matter 
term. Perturbing the neutrino emission in the transla¬ 
tional symmetry on the ring and in the emission direc¬ 
tions, we find the spontaneous breaking of these symme¬ 
tries in both normal and inverted mass hierarchies. As 
a consequence the flavor content and the lepton num¬ 
ber of the neutrino ensemble acquires seizable variations 
along different lines of sight. These findings are pre¬ 
sented as follows. In Sec. II we describe the features of 
our two-dimensional model. We discuss the equations of 
motion to characterize the two-dimensional flavor evolu¬ 
tion. We show how it is possible to solve this problem by 
Fourier transforming these equations, obtaining a tower 
of ordinary differential equations for the different Fourier 
modes. In Sec. Ill we present the numerical results of our 
study. We show how the breaking of the spatial symme¬ 
tries produce direction-dependent variations in the flavor 
content of the ensemble. Finally in Sec. IV we discuss 
about future developments and we conclude. 


II. TWO-DIMENSIONAL MODEL 
A. Equations of motion 

Characterizing the SN neutrino flavor dynamics 
amounts to follow the spatial evolution of the neutrino 


fluxes. For a stationary neutrino emission, the Equa¬ 
tions of Motion (EoMs) of the u space-dependent occu¬ 
pation numbers p(r, p) with momentum p at position r 
are 


V VrP = -f[Q, p] , (1) 

where we indicate with sans-serif vectors in flavor space, 
while for the ones in real space we use the bold-face. At 
the left-hand-side of Eq. © there is the Liouville opera¬ 
tor representing the drift term proportional to the neu¬ 
trino velocity v, due to particle free streaming. Note that 
we are neglecting external forces and an explicit time de¬ 
pendence of the occupation numbers. On the right-hand- 
side of Eq. o the matrix O is the full Hamiltonian that 
reads 




^ + oo 


Nj + 


dE'E' 


dw' 

(27r)‘ 


■q'{1 - V • v') 


where is the matrix of the mass-squared, responsible 
of the vacuum oscillations. The ordinary matter effects 
on neutrino flavor conversions is accounted by the matrix 
of charged lepton densities Ni. Finally, the neutrino- 
neutrino interaction potential is represented by the last 
term of the right-hand-side, where the integral in dV is 
on the unit sphere and the occupation numbers g' depend 
on r, E',\-'. Note that we use negative E to denote anti¬ 
neutrinos. 

In order to show the effect of spontaneous breaking 
of spatial symmetries we consider them to be emitted in 
a plane from a ring with radius r = R. We have then 
a two-dimensional model for which it is natural to use 
a system of polar coordinates to describe the neutrino 
position vector r = (r, (p) where r is the radius, (p € [0; 27r] 
is the polar angle, as shown in Fig. [TJ 

The neutrino velocity can be decomposed in the ra¬ 
dial {vr) and transverse {vt) component defined as v = 
(wri'Cib) = (cos0r,±sin0r), where 9r G [0,7r/2] is the 
angle between the radial direction and the one of the 
neutrino propagation (see, e.g., [s^), and the ± sign 
indicate a transverse velocity in the clock-wise (f+) or 
anti-clockwise (v-) direction with respect to the radial 
direction, respectively. We mention that in the recent 
multi-angle study |34| . where neutrinos emitted from a 
plane were considered, the range in the emission angles 
was 0 G [0,7r]. With this choice it is not necessary to 
distinguish clock-wise or anti-clockwise modes. However, 
in our work we preferred to use have the 9r G [0,7r/2] in 
order to start with a situation symmetric in the two emis¬ 
sion directions ± and show the effect of breaking of this 
discrete symmetry. Note that the local angle 6r would 
depend on the radius r. In order to avoid this effect, in 
the literature it is preferred to label the neutrino modes 
in terms of their emission angle G [0,7r/2] along the 
boundary at r = R. The two angles 9r and are related 
by di 


R sin 0R = r sin 9r ■ 


( 3 ) 
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FIG. 1: Two-dimensional model for the neutrino beams emit¬ 
ted from a ring with radius r = R. 


Furthermore, we introduce the angular variable u = 
sin0fl, u € [0,1]. With this choice the components of 
the neutrino velocity are [ 2 l| 


Vr = COs9r = \ll - 

R 

Vt = v± = ±sin^j^ = zb —u . 

r 


(4) 


We assume that the neutrino distributions at r = i? in 
the energy E and in the angular variables u and in the 
R, L directions can be factorized as 


F±{E, u) = F^E) X F^{u) X F± . (5) 


We assume the neutrino angular distributions to be flat 
in u and equal for all the flavors, i.e. Fy{u) = 1. 

Concerning the distributions in the clock-wise (-I-) or 
anti-clockwise (—) directions, we assume that these are 
given by 


F± 


(1 + /3t) 

2 + /?++/3_ ’ 


( 6 ) 


where the quantities ,5± <C 1 are introduced to slightly 
perturb the ± of a given u mode at the boundary. 

The neutrino number flux F^{E) at the ring is given 
by 


FAE) 


1 Li, 

^ttR^ {EA 


ME) , 


(7) 


where we have normalized the neutrino emission on a 
sphere with radius R. 


B. Two-flavor case 

In the following we will consider only a two-flavor sys¬ 
tem [veT^x) where x = fj,,T and we will describe the 


neutrino energy modes in terms of the two neutrino fre¬ 
quency uj = AnA I2 Eq, where Am? = m\ — m\ is the 
mass-squared difference. We have assumed a monochro¬ 
matic neutrino emission with E = Eq. In the two-flavor 
case the density matrices are projected over the Pauli 
matrices cr obtaining the polarization vectors in the usual 
way [T^, where we normalize the (anti)neutrino polariza¬ 
tion vectors to the difference of the anti-neutrino fluxes 
at the boundary. 

The Liouville operator on the left-hand-side of the 
EoMs [Eq. ([1])] assumes the form 


V • Vr 



v± d 
r d(j) 


so that the EoMs read (see also EH) 


( 8 ) 



v± d 
VrV d(j) 


p 


±,u 


— B + nA 

Vr 




(9) 


where we indicated with sans-serif the vectors in flavor 
space. In particidar, the unit vector B = 
points in the mass eigenstate direction in flavor space, 
such that B • 63 = — cos'd, where d is the vacuum mix¬ 
ing angle. For simplicity we neglect a possible matter 
effect, assuming that its only role would be to reduce 
the effective in-medium mixing angle, d <g; 1 j^. The 
neutrino-neutrino interaction terms has a multi-angle 
kernel (1 — Vp • Vq) which takes the form 


Vr 

R 


J d9'r[l — VrV'^ — vM]D' 


dd'n cos 


1 


- 1 - 


■ViVt 


VrV'j. 


VrV'j. 


D' (10) 


where D is the difference between the neutrino and anti¬ 
neutrino polarization vector of a given mode, and we used 
from Eq. ([3]) 


d9r 


R cos dfl 
r cos dr 




R cos dfl; , „ Rdu 

- ddR =-. 

r Vr r Vr 


( 11 ) 


In the large-distance limit (r ^ R) one can expand 
Eq. (flOl) obtaining 

J ■ (12) 


We note that for the case we are studying the self¬ 
interaction term declines as r“^, while in the SN case 
it declines as r~^. Considering the contribution of the 
clockwise (-I-) and anti-clockwise (—) modes in the pre¬ 
vious equation one gets 


vtM' = 

Vtv'M = T Uu’{D + ^u' - . (13) 
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Then, the neutrino self-interaction term in the large dis¬ 
tance limit r ^ R assumes the form 

f j /r/' 2 I /2,^ 

= fir du[[U +U )--- 

Jo ^ 

=F uu'{D+^u' - D-,«')] (14) 

where the =F refers to the ± modes respectively, and 


fir = [F,^iR) - F,^{R)] 


2r^ 


= 3.5 X 10^ km- 


15 MeV 
10^1 MeV/s 


_i /rY ( L 


rj y{E,^) 


10 km 
R 


(15) 


An equation analogous to Eq. (fTiHl can be written for the 
anti-neutrinos. 

One can define a conserved “lepton current” = 
(Lo,L) whose components are (see also 


Lo — 

F 1 

/ 2 (^+’"' 4- ^-,“') ■ ^ ’ 

(16) 

T — 

r 1 

/ du'Vr-{Pl + ,u' + P)-,u')-P> , 

Jo ^ 

(17) 

L* = 

F 1 

J du\vt\-iD^y, 

(18) 


where L is a two-dimensional vector (Lf.,Le), and 
• B ~ Dfl(L) „• From Eq. ([T]) one realizes that 
the lepton current satisfies a continuity equation 


cIqLq -t-Vj.'L — Vr'L — 0 , (m) 

where first equality follows since ^oLg = 0 having we 
assumed a stationary solution. Eq. (fT^ generalizes the 
lepton-number conservation law of the one dimensional 
case 0. 


where the initial flavor asymmetry is given by 

- Fr.. 


a = 


- Frr 


(23) 


The functions P±,«(r, 4>) are periodic in (j) with period 
27r. Their Fourier transform is then 


1 f 

P±,u,n{'^) = — J P±.«(T’»e- 


(24) 


so that 


p±.«(a<(>) 


+CO 

^ P±.«,n(r)e+™'^ 


n— — oo 


(25) 


The EoMs for the Fourier modes at large r ^ R as¬ 
sume the form 

d . R 

dr ’ ’ 

P ^ P±,ii,n 


-1-00 

4- 1*’’ F / d.u'[{u^ + u'^) 
„ Jo 


f2\ + ,u' ,n—j F ^ — ^u' 


j = -oc 

^UU yj_j)] X P 


We stress that it is enough to follow the evolution for pos¬ 
itive modes n > 0, since the P±,u(r, (j)) are real functions 
and therefore 


P±,u,r^ = P±,u,-n ■ (26) 

Once the evolution of the harmonic modes is obtained 
from Eq. (Ei), the polarization vector in configura¬ 
tion space can be obtained by inverse Fourier transform 
[Eq. EH)]. 


C. Equations of motion in Fourier space 

The differential operators in Eq. Q implies that the 
flavor evolution is characterized by a partial differential 
equation problem in r and <j)- In (see also [1^) it 

has been shown how it is possible to solve such a prob¬ 
lem by Fourier transforming the equations of motion with 
respect to the coordinate along which a perturbation is 
introduced. We assume a perturbation of the polariza¬ 
tion vectors at r = R with period 27r in (j) so that 

P±.u(I?,0) = P± + 2e,,5cos<^ , (20) 

where P^ is the unperturbed value of the polarization 
vector, and i5 <C 1 is the amplitude of the perturbation. 
Up to the small difference in the emission modes ± [see 
Eq. (IH] the initial values of the polarization vectors are 

Pliii) - (l + a)e, , (21) 

Pi(P) e, , (22) 


III. NUMERICAL EXAMPLES 

We present the results of the flavor evolution in the 
two-dimensional model described above. To calculate the 
li-ii interaction strength in Eq. m and the flavor asym¬ 
metry parameter in Eq. EH we use benchmark values 
often used in previous studies of self-induced neutrino 
oscillations (see, e.g. 0 ), i.e. we take as average ener¬ 
gies 

((E,J, (EpJ, (E^J) = (12,15,18) MeV , (27) 

while for the neutrino luminosities (in units of 10^^ erg/s) 
we assume 

L,, = 2.40 , Lp, = 2.0 , = 1.50 . (28) 

These values are typical of the early time SN accre¬ 
tion phase, and corresponds to an asymmetry parameter 
a = 1.34. As specified before we work in a single-energy 
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FIG. 2: Two-dimensional evolution of the 3-rd component P 3 
of the V polarization vector in the r-cf) plane, and its map on 
the bottom plane breaking only the ± symmetry. Upper plot 
refers to NH, while lower panel is for IH. 


FIG. 3: Two-dimensional evolution of the 3-rd component P 3 
of the v polarization vector in the r-<j) plane, and its map on 
the bottom plane breaking the azimuthal invariance for the 
± modes and the translational symmetry on the ring. Upper 
plot refers to NH, while lower panel is for IH. 


scheme, where we take as representative vacuum oscilla¬ 
tion frequency the one corresponding to the average of 
the V ensemble with the emissivity parameters chosen 
above (see [3|)- Namely we take lo = 0.68 km“^. Con¬ 
cerning the neutrino oscillation parameters we choose a 
small mixing angle d = 10“^. Moreover, we assume 
Nu = 100 modes for the angular variable u € [0; 1] in 
order to have numerical convergence of the results and 
to avoid spurious instabilities due to few angular modes. 

It is known that forcing the ± symmetry (taking 
/3+ = /3_ = 0 in Eq. ®) and the translational sym¬ 
metry on the ring (taking (5 = 0 in Eq. (l20ll i the ensem¬ 
ble is stable in normal mass hierarchy (NH, Am^ > 0) 
while in inverted mass hierarchy (IH, Am^ < 0) it ex¬ 
hibits large bimodal flavor changes in the form of pair 
conversions v^Vf. VxVx Q- If perturb the ± sym¬ 
metry taking small seeds /3+ = —/3- in the distributions 
of Eq. (|ni) the system now exhibits the analogous of the 
so-called multi-azimuthal-angle (MAA) instability of the 
bulb model [2l|, [23,113 . The result of the flavor evolution 
is shown in Fig. [2] where it is represented the behavior 
of the 3-rd component of the (integrated over u) anti¬ 
neutrino polarization vector P3 = 1/2(P+ -|- P_) in the 
(r, (j)) plane for NH (upper panel) and IH (lower panel). 
In these numerical examples we have chosen /?+ = 10“^. 
The most striking effect of the MAA instability is that 


now also NH exhibits flavor conversions at r ^ 60 km. 
The choice of the initial seed /3± determines the onset 
radius of the flavor conversions: the largest the seed, 
the earliest flavor conversions start. In IH (lower panel) 
flavor conversions start as in the ± symmetric case at 
r > 50 km and MAA has a minor impact on the flavor 
evolution. From the Figure we realize that the behav¬ 
ior of the flavor conversions is uniform in the (j) variable 
since the translational symmetry has remained unbroken. 
Indeed we have solved only the EoMs [Eq. (l26ll ] for the 
n = 0 Fourier mode in this case. 

The next step is to perturb also the translational sym¬ 
metry on the ring assuming a seed 5 in the longitutudinal 
distribution of the polarization vectors on the boundary 
[see Eq. (l20l) ]. In this case we consider the evoltion of the 
first N=I00 Fourier modes in Eq. (l26)l . In this way we 
are sensitive to variations occuring at an angular scale 
A(j) > 3°. Results are shown in Fig. [3] with the same for¬ 
mat of the previous Figure. We used as seed to break the 
^-symmetry d = 3 x 10“^. In both NH (upper panel) and 
IH case (lower panel) flavor conversions start as in the 
translational invariant case, i.e. the planes of common 
oscillation phase are flat in (j) direction. However this be¬ 
havior is not stable. In the NH case around r ~ 100 km 
something occurs: The P3 component is no longer flat in 






6 



X (km) 


FIG. 4: Component Lr of the vector lepton number L in carte¬ 
sian coordinates in NH (upper panel) and IH (lower panel), 
respectively. 
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FIG. 5: Contour plots of the first 100 Fourier modes |Pn| (in 
logarithmic scale) in the plane n-r in NH (upper panel) and 
IH (lower panel) respectively. 


(p, while it starts to acquire notable variations (~ 20 %) 
at different longitude. In IH after flavor conversions de¬ 
velop, also the translational symmetry is perturbed at 
r ^ 100 km, with variations at different (p with values up 
to 30 %. 

In Fig. [4] we represent the component Lr of the vector 
lepton number L [Eq. dlSp ] in cartesian coordinates 

X = rcos(p , 

y = rsinc/) . (29) 

We realize that when the spherical symmetry is broken, 
the lepton number acquires significant variations in dif¬ 
ferent directions at a given r with respect to the initial 
uniform value = a = 1.34. In particular one finds 
~ 20 % variations. 

In order to clarify better this flavor dynamics, in Fig. [5] 
we show a contour plot representing the growth of the 
different Fourier modes | P„ | (in logarithmic scale) in the 
plane n-r for NH (upper panel) and IH (lower panel). We 
consider the evolution of the first N = 100 modes. We 
realize that the breaking of the translational symmetry 
at r ~ 60 km in NH corresponds to the rapid excitation 
of the n > 0 harmonics that reach values |P„| < 10“^. In¬ 
stead in IH for r > 50 km the modes start to get excited 


and can grow to |P„| ^ 10“^'® at r > 150 km. This ex¬ 
plains why the effect of the breaking of the translational 
symmetry is more pronounced in IH rather than in NH. 


IV. CONCLUSIONS 

We have considered a simple two-dimensional toy- 
model, namely neutrino beams emitted from a ring, to 
point-out the effect of spontaneous breaking of axial and 
spherical symmetries in the self-induced flavor conver¬ 
sions of SN neutrinos. 

We found that if the slightly perturb the space symme¬ 
tries on the boundary, these perturbation seeds are dra¬ 
matically amplified altering the flavor conversions found 
in a symmetric model. Therefore, the flavor content of 
the self-interacting SN neutrinos would acquire signifi¬ 
cant direction-dependent variations. These results are 
qualitatively similar to what we found in the planar 
model we studied in H- Our findings suggest that the 
characterization on the flavor conversions obtained before 
should be critically reconsidered, including these sponta¬ 
neous symmetry breaking effects. In order to have a re¬ 
alistic characterization of the possible SN neutrino spec¬ 
tra our simple toy model should be improved on differ¬ 
ent aspects. In particular, one should extend this model 
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to a realistic three-dimensional spherical case. In this 
case one would have the possibility to break the spherical 
symmetry in both longitudinal and latitudinal directions. 
Moreover, in order to get our numerical solution we have 
considered N = 100 Fourier modes. In this cases we have 
not found the presence of flavor converions at lower radii 
than in the spherically symmetric case. However, in 
it has been shown that harmonics with sufficiently high 
n could become unstable also at low-radii. Increasing N 
to 500 we have not found any sizeble change in the onset 
of the flavor changes and in the subsequent flavor evolu¬ 
tion in the non-linear regime. However, it remains to be 
seen if with a much higher number of harmonics low-radii 
effects could occur. In this regard, it would be useful a 
stability analysis performed along the lines of [s^- 
Continuous energy spectra should also be taken into 
account to understand how the spectral splitting fea¬ 
tures found in the bulb model would be modified in 
this case. The role of matter effects that would sup¬ 
press self-induced flavor conversions during the accre¬ 
tion phase should also be investigated. The final goal 
would be to study of the self-induced neutrino flavor con¬ 


versions in realistic multi-dimensional supernova models 
accounting for largely aspherical neutrino emission and 
matter profiles. This objective is particularly timely now 
since in the last recent years, SN model simulations have 
experienced several breakthroughs. After ID [s^ and 
2D m odels, the forefront has reached 3D SN simula¬ 
tions Therefore, it seems the perfect juncture to 

connect realistic SN simulations with nonlinear neutrino 
oscillations. This open issue makes compulsory the need 
for further dedicated studies to fully clarify the fascinat¬ 
ing behavior of the interacting neutrino field. 
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